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Abstract 

Following the discovery of the cosmic microwave background, the 
hot big-bang model has become the standard cosmological model. In 
this theory, small primordial fluctuations are subsequently amplified by 
gravity to form the large-scale structure seen today. Different theories 
for unified models of particle physics, lead to different predictions for 
the statistical properties of the primordial fluctuations, that can be di- 
vided in two classes: gaussian and non-gaussian. Convincing evidence 
against or for gaussian initial conditions would rule out many scenarios 
and point us towards a physical theory for the origin of structures. 

The statistical distribution of cosmological perturbations, as we 
observe them, can deviate from the gaussian distribution in several 
different ways. Even if perturbations start off gaussian, non-linear 
gravitational evolution can introduce non-gaussian features. Addition- 
ally, our knowledge of the Universe comes principally from the study of 
luminous material such as galaxies, but galaxies might not be faithful 
tracers of the underlying mass distribution. The relationship between 
fluctuations in the mass and in the galaxies distribution (bias), is often 
assumed to be local, but could well be non-linear. Moreover, galaxy 
catalogues use the redshift as third spatial coordinate: the resulting 
redshift-space map of the galaxy distribution is non-linearly distorted 
by peculiar velocities. Non- linear gravitational evolution, biasing, and 
redshift-space distortion introduce non-gaussianity, even in an initially 
gaussian fluctuation field. 

I will investigate the statistical tools that allow us, in principle, to 
disentangle the above different effects, and the observational datasets 
we require to do so in practice. 
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1 Introduction 



Until recently in cosmology, non-gaussianity has been a synonymous of non- 
linearity; but, in the last 5 years or so, more and more objects like the 
galaxy of [1] at redshift 5.6 have been found. For the first time a galaxy has 
been found at higher redshift than the most distant known quasar. More 
recently, a galaxy at redshift almost 7 has been found [2]. The standard 
"inflationary" cosmological model with gaussian initial conditions predicts 
that these objects should be very rare. It is becoming increasingly difficult 
to accommodate the existence of so many high-redshift galaxies under the 
assumption that non-gaussianity is equivalent to non-linearity, that is pos- 
tulating gaussian initial conditions. Non-gaussianity does not necessarily 
imply non-linearity: there might be some primordial non-gaussianity and it 
is necessary to "find a way" to distinguish the two effects. 

2 Non-gaussianity due to non-linearities 

Let us start by assuming gaussian initial conditions and investigate the ef- 
fects of non-linearities. We define the fractional density contrast 5 as 5p/p, 
where p is the mean density. The probability distribution of 6 starts off 
symmetric around zero, with negligible tails for \5\ > 1. Non-linear gravi- 
tational evolution skews the distribution towards high densities: this is due 
to the fact that underdense regions cannot become more empty than the 
void (5 > —1) while overdense regions can accrete matter arbitrarily (no 
upper limit on 6). This is not the only process that can skew an initially 
gaussian distribution. The mass in the Universe is mainly dark matter and 
cannot be observed directly: only galaxies can be observed, but mass and 
galaxy distributions may not be identical: the idea that galaxies are biased 
tracers of the mass distribution was introduced in the early eighties [3]^ and 
has featured strongly in large scale structure (LSS) studies. In general, bias 
must alter the statistics of any underlying matter distribution, otherwise 
6 < —1 for the galaxy field, which corresponds to a negative galaxy density. 
In different bias schemes suggested in the literature, the relation between 
the galaxy and the mass fluctuation fields {6g and 6 respectively) has been 
taken to be local, non-local, eulerian, lagrangian, stochastic etc... 

In what follows we will assume that Sg{x) = F[6{x)], that is the bias, is 

^Although the fact that galaxies of different morphologies have different spatial dis- 
tributions and they cannot all be good tracers of the underlying mass distribution, was 
known much before the introduction of the concept of bias (e.g. [4]). 
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a local eulerian function of the underlying mass field. Furthermore we will 
assume (following [5]) that F can be expanded in Taylor series and we will 
truncate the expansion to the quadratic term: 



60 is unimportant and simply ensures that {6g) = 0. This non-linear opera- 
tion on the matter field introduces some skewness, i.e. some non-gaussianity. 
As first suggested by Fry [6] it is possible to disentangle the two non-gaussian 
contributions, non-linear gravity and bias, by looking at higher-order cor- 
relations in the mildly non-linear regime. In particular, if the initial fluc- 
tuations are gaussian and cosmological structures grow by gravitational in- 
stability, the three-point correlation function is intrinsically a second-order 
quantity^ and is detectable in the mildly non-linear regime. If then bias 
can be expressed as in equation (|l]), it is possible to show that a likelihood 



where 6k is the Fourier transform of 5(x). Due to the presence of the Dirac 
delta function 6^, the bispectrum can be non zero only when the three k 
form a triangle. 

In practice, the higher-order statistic (the bispectrum) exploits the fact 
that gravitational instability skews the density field as it evolves, creating 
sheets and filament-like structures reminiscent of the Zeldovich pancakes. 
Non-linear bias also introduces skewness but does so by shifting the iso- 
density contours up and down, without modifying the shape of the struc- 
tures. These two effects can be disentangled by using different triangle 
shapes for the bispectrum. 

There are several advantages in performing this sort of analysis in Fourier 
space, most of them are the same advantages of the power spectrum over the 
two-point correlation function. We will recall here only the following: the 
estimates of power on different scales can be made uncorrelated; it is easy to 
deal with the error estimate; and more importantly it is easy to distinguish 

^In the quantity S, assumed to be small. 

^The likelihood method can easily be generalized to measure the lagrangian [7] and 
stochastic (e.g. [8]) bias parameters. 
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between linear, mildly non-linear and highly non-linear scales. In the mildly 
non-linear regime the bispectrum is given by: 



) = {27rfPg{h)Pg{k2)[ciJ{ki,k2)+C2] 5^(ki+ k2+ ks)+cyc. 

(3) 

where Pg denotes the galaxy power spectrum, J(ki, is a known function 
of the two A;-vectors and 

ci = ^ C2 = I . (4) 

In the absence of bias (ci = 1, C2 = 0) it would be easy to isolate the non- 
gaussianity generated by gravitational instability. Nevertheless, even in the 
presence of bias, it is possible to disentangle gravity from bias. Equation 
(^) is in a form suitable for a likelihood analysis for the two bias parameters 
via ci and C2 [9], once the covariance is knownQ. A generating functional 
approach to calculate analytically the N-point function and therefore the 
covariance for the bispectrum was introduced in [10,11]. The performance 
of the method has been tested on biased and unbiased N-body simulations^, 
with very promising results^ [11] for the application to forthcoming galaxy 
redshift surveys such as SDSS and 2dF. An estimation for the expected error 
achievable from present galaxy surveys, yields an error on ci of about 100% 
[11] and is therefore not particularly useful. 



3 Real world issues 

Of course, reality is always more complicated: in a realistic galaxy survey 
several complications arise due to the presence of shot noise and selection 
function, but more importantly due to redshift space distortions. Galaxy 
surveys in fact use the redshift as the third spatial coordinate. The redshift 
would be an accurate distance indicator in a perfectly homogeneous Uni- 
verse; but the universe is clumpy, inhomogeneities perturb the Hubble flow 
and introduce peculiar velocities. The resulting redshift-space map of the 
galaxy distribution is thus distorted along the line of sight, and the nature of 
this distortion is intrinsically non-linear. On large scales the coherent inflow 

''The bispectrum is a three point quantity, its covariance is a six-point quantity 
(pentaspectrum) . 

^The N-body simulation was provided by the Hydra-consortium and produced using 
the code [12]. 

®For a different approach see J. Frieman contribution in this volume. 
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into overdense regions introduces a squashing effect in the redshift map (the 
great wall effect), on smaller scales, the virialized highly non- linear struc- 
tures appear elongated along the line of sight (the fingers-of-God), heavily 
contaminating the mildly non-linear regime were most of the signal for the 
bispectrum comes from. In [13] we showed that, with an accurate model- 
ing of redshift space distortions in the distant observer approximation, it 
is possible to disentangle the effects of redshift-space anisotropics from the 
bias and gravitational effects. This is achieved by combining a second-order 
perturbation theory description of the coherent inflow (e.g. [14]) with an ex- 
ponential velocity dispersion model, and discarding the A;-modes where the 
contamination from highly non-linear structures is too big to be successfully 
modeled. An alternative approach hss also been explored (e.g. [15]). 
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Figure 1: Likelihood for ci = 1/62 and 02 = b2/h\, (where bi andb2 are linear 
and quadratic bias parameters) for an unbiased simulation (bi = 1, b2 = 0) 
in redshift space. The solid and dotted lines show the 1 and 3 a confidence 
levels respectively. This fi,gure shows that it is possible to disentangle the 
non-linear gravitational instability, bias and redshift space distortion signals, 
and also that from future galaxy redshift surveys such as SDSS and 2dF, the 
bias could be known with an accuracy better that 10%. 



The result of the likelihood analysis performed on an redshift space un- 
biased simulation of lOO/i^^ Mpc side is shown in figure 1. 

It shows not only that it is possible to disentangle the non-linear gravi- 
tational instability, bias and redshift space distortion signals, but also that 
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from future galaxy redshift surveys such as Sloan digital sky survey (SDSS) 
and the Anglo- Australian two-degree field (2dF), the bias could be known 
with an accuracy better that 10%. Before achieving this goal however, there 
are several other issues to deal with, that arc potentially serious for any 
Fourier based technique. In particular, these are the mask, the difficulty of 
obtaining redshifts for close pairs of galaxies, the holes arising from bright- 
star drills, and the variable completeness. We have investigated and quan- 
tified these effects on a simulated catalogue and concluded that a 10% error 
on the bias parameter could be achieved from the 2dF survey [16]. 

3.1 Bypassing the redshift-space distortions 

Two-dimensional surveys avoid the redshift-space distortion problems, but, 
in principle, contain less information. However, because of the smaller obser- 
vational effort required, these can contain a much larger number of objects. 
For example the APM survey at present contains 10^ galaxies, the DPOSS 
catalogue will have 50 million galaxies and the SDSS will provide us with a 
two-dimensional map of 10^ galaxies. 

To treat the projection of higher-order correlations on the celestial sphere, 
the spherical nature of the distribution cannot be ignored. Spherical har- 
monics are eigenfunctions for the two-dimensional surface of the sphere and 
therefore are the natural basis describing a two-dimensional random field on 
the sky. When comparing with Fourier space analysis we have that 

4 af (5) 

and in particular for the bispectrum 

<5^(kik2k3)^(^4t.'^3) (6) 

where on the RHS we have the three-J symbol. To perform the same sort of 
analysis as the one illustrated in sections 2 and 3 an expression that relates 
the 3D bispectrum to the projected one in spherical harmonics is needed: 



116 /(24 + I) (2^2 + 1) (2 4 + 1) 
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d/cidA;2Z^i+^^A;?A;i*^,(A;i)^'^,(A;2) ^ ^^«+^^(-l)^S^(A;l, A;2)(24+l)(24+l)^x 
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/ drr-i.{r)n,{k,r)u,{k,r) (^,^ ^) (^i ^l) {Wn} + eye] 

This expression is quite eomplicated: the derivation and the detailed expla- 
nation of it can be found in [17]. For the purpose of this contribution we 
only have to notice that it is an exact expression relating the spherical har- 
monics projected bispectrum (a™^a^^a™^) to the 3D bispectrum expressed 
through its Legendre coefficients i?£(kjkj). (fcj) is a known function of 
the selection function, denotes the spherical Bessel function and {. . .} 
denotes the Wigner 6-J symbol. 
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Figure 2: Likelihood contours for degenerate triplets configurations. The 

two levels are the 1—a and 3-(t confidence levels and the + indicates where 
the true value for the parameters lies. Perturbation theory breaks down at 
i bQ. Adding other configurations does not help much: this is about the 
best result obtainable from projected surveys. 

By performing a likelihood analysis to measure the bias parameter on a 
all sky simulation with the APM selection function, we find that the results 
are not encouraging for projected catalogues (see figure 2). 

It is therefore preferable to undertake the bispectrum study of 3D galaxy 
redshift surveys such as SDSS and 2dF using the method described above. 
The good news is that the exact expression for the projected bispectrum 
in spherical harmonics has applications in a variety of areas such as cosmic 
microwave background (CMB) and gravitational lensing studies. 
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4 Primordial non-gaussianity: CMB vs. LSS 



Up to now we have assumed gaussian initial conditions. However, among 
the theories for structure formation, only inflation predicts a nearly gaussian 
distribution for the primordial fluctuations, with deviations from gaussianity 
which are calculable, small and dependent on the specific inflationary model 
(e.g. [18,19,20]). In other models such as non-standard inflation or topo- 
logical defects models initial conditions are non-gaussian. CMB and LSS 
data will shortly improve dramatically: it is therefore timely to ask which 
of the CMB or LSS will provide a better probe of the nature of primordial 
fluctuations. The advantage of looking at CMB maps is that the fluctuation 
distribution should be close to the primeval form, but the disadvantage is 
that the amplitude of fluctuations is small and there are foregrounds and 
other effects to account for. The advantages of looking at LSS is that the 
signal has been amplified by gravity. This is however also a disadvantage 
because gravity skews the distribution. Non-linear gravity, bias and redshift 
space distortions might completely swamp the primordial signal. Following 
[21], as a discriminating statistic we will use the bispectrum, but we will 
start by considering the skewness as an example to illustrate some of the 
effects. 

The skewness is defined as 

For a gaussian field the skewness is zero, while for an initially gaussian field 
evolved under gravitational instability to second order in 6, the skewness 
becomes 34/7 and is constant in timeQ In what follows for CMB related 
calculations, we will assume an Einstein de Sitter Universe, this assump- 
tion is justified because we shall be concerned with factors of 10 while the 
cosmology can change the results only by factors of order unity. Suppose 
the initial conditions are very close to gaussian, but with a small primordial 
skewness, parameterized by S^^z = 1100) at recombination. The effect on 
LSS win be (e.g. [22]): 

53 = ^3,0 + 34/7 + .^ (9) 
where J- includes a complicated dependence on the three- and four-point 



^The value 34/7 is strictly true only for an Einstein de Sitter Universe, but the skewness 
does not depend strongly on cosmological parameters. 
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function that we will ignore for the moment and 53^0 scales as: 

the primordial skewness redshifts away. We can then make a thought ex- 
periment: assume we know the real space position of every particle in the 
whole Hubble volume. The smallest error for the skewness, that is the 
smallest 5*3^0 detectable, on 20 Mpc scales is 5*3^0 ~ 10~^, which implies 
5*3 (z = 1100) ~ 10. We can repeat the exercise for the CMB where now, 
for consistency, we consider the smallest detectable skewness on 0.2° scales, 
obtaining 53(2; = 1100) ~ few < 10. 

This example already shows that CMB seems to be more sensitive to pri- 
mordial deviation from gaussianity than LSS, but we will now proceed more 
accurately by considering the bispectrum: in fact the bispectrum contains 
more information than the skewness and has all the advantages of being 
a Fourier space quantity. In the absence of bias, the LSS bispectrum in 
second-order perturbation theory for non-gaussian initial condition is: 

5(ki,k2,k3) = So(ki,k2,k3) + 2J(ki,k2)P(A:i)P(fc2) + eye. 

+ j d3A:J(k',k3-k')r^(k',k3-k',ki,k2)+eye. (11) 

Bq is the primordial bispectrum evolved linearly and corresponds to 6*3^0 
of equation ^ the second term is the usual gravitational instability bispec- 
trum and correspond to the 34/7 term in equation ^ T'^ denotes the Fourier 
counterpart of the connected four point function and the integral term cor- 
respond to T of equation ^. We then parameterize the LSS bispectrum 
as: 

B = P(A:i)P(/c2)[2J(ki, k2)ei + e2] + eye. (12) 

because we know how to estimate eiand e2 from LSS studies (section 2). In 
the very idealized case where the real space position of every particle in the 
SDSS volume was known, the minimum ei and e2 detectable would be re- 
spectively ei ~ 10"'^ and e2 ~ 10~^ ignoring all the real world complications 
of shot noise, selection function etc... 

On the CMB side, the bispectrum for realistic non-gaussian models is 
given by (e.g. [23]): 

Bt,i,i, = f{h,e2j3)a[Ce,Ce,+cyc.] (13) 

where a is the amplitude, / is a known function of the i s and Ci denotes 
the CMB power spectrum. The minimum error on the amplitude a can 
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be obtained via tlic Fisher information matrix: 




E 




(14) 



m\m2mz 



N{mi,£i) 



where n = 1/2 and N{mi,ii) is the number of non-zero terms Hke Ci-^Ci^Ci.,^ 
in the covariance and ranges from 1 to 30. Here we neglect partial sky 
coverage effects, and by constraining i < 100 pixel noise and small angular 
scale effects are negligible. 

In [21] we investigated the LSS and CMB bispectrum as a discriminating 
statistic for several physically motivated non-gaussian models. There is an 
infinitude of deviations from gaussianity and one cannot address them all: 
we thus restrict ourselves to physically motivated models where the non- 
gaussianity can be dialed from zero (the gaussian limit) and is assumed to 
be small. In particular we consider the non gaussianity parameterized by: 



where (f) denotes a gaussian field and for the moment we will assume that 
$ is the gravitational potential; the non-gaussianiy parameter is e, that is 
zero for a gaussian field, ~ 1 for standard inflation, but can be as big as 
^ 20 for some non-standard inflationary models. The CMB effect is given 
by 2€/Asw = o: where Asw is the Sachs- Wolfe coefficient 1/3. 

It is possible to see [21] that, if the CMB distribution from the future 
satellite missions turns out to be consistent with gaussian, the smallest e 
allowed would be ^ 20. 

The LSS effect is C2 = 62/^1 + 10^^e and the T'^ contribution is negligible. 
By substituting the minimum e measurable from CMB in this expression 
and ignoring bias, we obtain C2 = 10~^: about two orders of magnitude 
smaller that the minimum C2 detectable from LSS even in the most idealized 
conditions. 

In [21] also other physically motivated non-gaussian models have been 
considered, the result is qualitatively always the same: if future CMB maps 
are consistent with the gaussian hypothesis then any non-gaussianity seen in 
the LSS bispectrum is due to non-linear gravity or bias, and we know how 
to disentangle the two. 

5 Looking at smaller scales 

In practice, CMB studies can be affected by noise and foreground and, 
more importantly, there might be models in which non-gaussianity is present 



^ = </' + e(0'-(0')) 



(15) 
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mainly on LSS or galaxy scales, which are not fully accessible with CMB ex- 
periments. I will therefore investigate another two ways to detect primordial 
non-gaussianity on scales smaller than CMB ones. 

5.1 Detecting non-gaussian initial conditions from large-scale 
structure 

It is possible to bypass the contamination due to non-linear clustering and 
discriminate between gaussian and non-gaussian initial conditions by using 
higher-order statistics in LSS studies such as the trispectrum -that is the 
connected four-point correlation function in Fourier space [24]. This quan- 
tity has the advantage of having a rather simple growth rate, with no com- 
plicating contributions from non-linear gravity in second-order perturbation 
theory and that the analysis depends on cosmology and bias only through 
the measurable quantity^ /?. The departures from gaussian statistics can 
be parameterized by introducing the quantity H, which is the fractional ex- 
cess of the 4-point function over the gaussian (disconnected) trispectrum. 
This quantity, in specific cases, can give us a meaningful measure of 'non- 
gaussianity': for mildly non-gaussian fields, a measurement of H can reliably 
be made. For highly non-gaussian fields, the gaussian hypothesis can be re- 
jected, but the measurement of H will be unreliable. Following [24] it is 
possible to deal with redshift-space distortions, biasing, spatially varying 
selection function and shot-noise. Figure 3 shows the minimum analysis 
for the parameter H from a redshift-space unbiased CDM-like simulation. 
By applying this method to future galaxy surveys, such as the SDSS, it will 
be possible to place tight constraints on initial departures from gaussian 
behavior. 

5.2 The abundance of high-redshift objects as a probe of 
non-gaussian initial conditions 

LSS probes scales much larger than galaxies but smaller than those acces- 
sible by CMB observations, and probes the present-day Universe at z = 0; 
conversely CMB maps probe the Universe at redshift z 1100 and even 
larger scales. The abundance of cosmological structures at redshift in be- 
tween these two ends and in particular at z > 1, contains also vital infor- 
mation about the nature of primordial fluctuations due to the fact that one 

®The quantity /3 arises naturally when studying the large-scale squashing effect of 
structures in redshift-space maps. It is defined as /3 ~ Qq'^ /b where b is the linear bias 
parameter. 
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Figure 3: Minimum analysis for the parameter H from a redshift-space 
unbiased CDM-like simulation. The analysis is largely bias independent. 
The true value for H is 0, and it is nicely within the 1-a level (minimum 
+0.5.). The quantity H , in specific cases, can give us a meaningful measure 
of 'non-gaussianity ': for mildly non-gaussian fields. 

is probing the tail of the distribution. Figure 4 shows that the effect of 
a small non-gaussianity is dramatic on the tails: high peaks are amplified 
much more than low ones and deep through can become local maxima. 

To extract this information, the Press-Schechter (PS) formalism [25] 
needs to be extended to non-gaussian initial conditions. The PS is an an- 
alytical model to calculate the mass function (that is the number of object 
per unit mass at a given redshift per unit volume), within an appropriate 
theoretical model (see R. Sheth and S. Shandarin contributions in this vol- 
ume). The key ingredient is the probability density function (PDF) for the 
smoothed dark matter field V(5m)] in fact, the number density of objects 
above a given mass M (corresponding to a smoothing radius R) at a given 
redshift (the mass function)\s proportional to the quantity P>j^: 

P>sS^m) = / V{SM)dSM (16) 

JSc 

where 5c is the threshold overdensity for the object to collapse, P^sd^M) is 
evidently a function of the redshift of formation (or collapse) of the object Zc, 

this redshift dependence is enclosed in 5c- 5c.{z) = Ac/D{z). Here D{z) is the 
linear growth factor, which in turn depends on the background cosmology, 
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Figure 4: The effect of a small non-gaussianity are dramatic on the tail of the 
distribution. Assume 4> here denotes a one dimensional gaussian fluctuation 
filed (black line); the dashed line is given by cj) + (j)'^ . High peaks are amplified 
much more than low ones and deep through can become local maxima. 



and Ac is the linear extrapolation of the overdensity for spherical collapse; 
Ac is traditionally taken to be ~ 1.68 but other values have also been used 
(see S. Shandarin and R. Sheth contributions). 

For gaussian fields, V{5m) is of course well known, but needs to be 
computed for non-gaussian initial condition. 

Given a physically motivated parameterization of primordial non gaus- 
sianity, we set off to calculate the PDF for the smoothed dark matter field 
analytically. Of course one could evaluate the PDF from numerical simula- 
tions, but this approach is plagued by the difficulty of properly accounting 
for the non-linear way in which resolution and finite box-size effects propa- 
gate into the statistical properties of the non-gaussian field |^. 

As before we parameterize the non-gaussianity as in equation 15 where 
<1> can be the potential or the density field. To properly deal with the 

^For example imagine computing the power spectrum of a non-gaussian field ip ~ 
cj> + (f)^ ~ {(fy^) where (j) is gaussian with a power-law power spectrum P^. The power 
spectrum of 4' involves computing the convolution of two P^. This convolution is an 
integral over k from to infinity. When performing the operation on a simulation box the 
final result would be as if the integral was truncated at > 2ti/L and k < 2tt/1 where L 
is the side of the box and I is the grid resolution. 
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smoothing we use a path-integral approach in the calculation of the PDF: 



cf>R{x) + e / cISFrUx - y|)</)2(y) -C- 6r{x) 



(17) 



where the R subscript denotes the smoothed quantity, 6^[. . .] the Dirac delta 
function, (. . .) the ensemble average and F is defined through its Fourier 
transform (FR(k) = WR{k)T{k)g{k) with Wr the smoothing, T the transfer 
function and g = 1 for the density ov g = —2/3{k/Ho)^Q^l^ for the po- 
tential) . In (f)R smoothing and transfer function are easily accounted for, 
but in the non-gaussian part, the presence of the smoothing in Fr{\x — y\), 
makes the quantity non-local. The Dirac delta function can be expressed in 
its integral representation and the ensemble average can be written as an 
integral over all (p configurations weighted by the gaussian probability den- 
sity functional. In this way we can express an unknown quantity in terms 
of all known quantities and we are able to obtain the non-gaussian PDF for 
the smoothed field analytically. The details of the derivation can be found 
in [26]. The main result is that, for mildly non-gaussian initial conditions 
with small positive skewness S^^r, the threshold for collapse 6c is lowered, 
in particular [26]: 



6c{Zc) > 6c{Zc) 



1 5-<5c(2;c) 



(18) 



where S-^^r is directly proportional to the non-gaussianity parameter e of 



equation 15. By lowering the threshold for collapse rare objects will form 
more easily, and this has a huge impact on the tails of the distribution as 
shown in figures 5 and 6. 

In particular figure ^ shows that the non-gaussianity of equation 15 ap- 
plied to the density field with e ~ 10~^, can change the number density 
of objects of mass M ^ 10"*^^ Mq that collapse at redshift Zc = 8, by two 
orders of magnitude. Conversely, it is clear from figure |^ that non-gaussian 
mapping |l^, applied to the potential field, has dramatic effects on cluster 
scales. Observations of clusters with Zc > 2 and M > 10^^ could put 
some constraints on inflationary models [27]. 



5.2.1 A worked example 

Up to date 6 galaxies with confirmed spectroscopic redshifts have been ob- 
served with redshifts 5 < z < 7. The observed comoving density N is for a 
Qq^^ = 0.3 , A = 0.7 (ACDM) Universe is iV > 8.3x 10-^(/i-iMpc)-3,[7V > 
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Figure 5: Effects of the non-gaussianity of equation^^ applied on the density 
fluctuation field. The figure shows only galaxy mass scales: the solid line is 
the mass function for gaussian initial conditions, all the other lines have 
non-gaussianity parameter 10~^ ^ e ^ 10^^. For Zc = 8, a small non- 
gaussianity parameter e ~ 10~^ increases the mass function by two orders 
of magnitude at about M ~ 10^^ Mq. 



3.6 X 10~^(/i~^Mpc)~^ for an Einstein-de Sitter Unverse]. Their masses are 
very uncertain, but some estimate can be obtained with simple arguments 
about their observed star formation rate, these estimates can then be com- 
pared with Lyq, width observations [26] . The gaussian ACDM model predicts 
> 5.2 X 10^^ (h^^Mpc)"^ , a factor ~ 20 fewer objects, while the Einstein- 
de Sitter model predicts > 10~^(/i~^Mpc)~''^: a factor 10^ fewer objects! 
Only e ~ 10~^ in the density is needed to reconcile ACDM predictions with 
observations. Alternatively, this discrepancy could be explained postulating 
an error of about a factor 4 in the mass determination. As larger telescopes 
such as NGST get on line it will be possible to determine masses more ac- 
curately and thus constrain the amount of primordial non-gaussianity on 
galaxy scales. 



6 Conclusions 

We have shown that non-gaussianity does not necessarily mean non-linearity, 
but it is possible to distinguish different kinds of non-linearity: e.g. bias, 
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Figure 6: The non gaussianity of equation It applied to the gravitational 



potential field has dramatic effects on cluster scales. In particular observa- 
tions of clusters with ~ 2 — 3 and M > W^^Mq could put constraints on 
inflationary models. The absolute value for e here is relatively big, but the 
deviation from gaussianity is still small: the skewness here is of the same 
order of magnitude as in figure 5. 



gravitational evolution, redshift space distortions. For the physically moti- 
vated non-gaussian models we considered, it turns out that CMB bispectrum 
is better that LSS bispectrum to detect primordial non-gaussianity: if the fu- 
ture CMB missions will produce maps that are consistent with the gaussian 
hypothesis, any non-gaussianity seen in the LSS bispectrum can be unam- 
biguously attributed to the effects of non-linearities. Thus, if this is the case, 
from on-going LSS surveys such as SDSS and 2dF we will be able to know 
the bias with few % accuracy. We have also shown that, to measure the bias 
parameter, ongoing 3D surveys are much better that 2D ones, even with 
full sky coverage, but the method developed has applications in different 
areas such as CMB and gravitational lensing studies. To conclude, we have 
seen different ways to disentangle primordial non-gaussianity from effects 
of non-linearity: CMB bispectrum, LSS trispectrum and the abundance of 
high-redshift objects such as galaxies and clusters. These methods probe 
the Universe a different scales and at different times and in addition to that 
they are sensitive to different moments of the distribution. We should there- 
fore conclude that these methods are complementary rather than mutually 
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exclusive. 
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